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Abstract. Spitzer/IRAC images of extended emission provide a new insight on the nature of small dust particles in the Galactic 
diffuse interstellar medium. We measure IRAC colors of extended emission in several fields covering a range of Galactic lati- 
tudes and longitudes outside of star forming regions. We determine the nature of the Galactic diffuse emission in Spitzer/IRAC 
images by combining them with spectroscopic data. We show that PAH features make the emission in the IRAC 5.8 and 8.0 /jm 
channels, whereas the 3.3 yum feature represents only 20 to 50% of the IRAC 3.6 fim channel. A NIR continuum is necessary 
to account for IRAC 4.5 yum emission and the remaining fraction of the IRAC 3.6 yum emission. This continuum cannot be 
accounted by scattered light. It represents 9% of the total power absorbed by PAHs and 120% of the interstellar UV photon 
flux. The 3.3 yum feature is observed to vary from field-to-field with respect to the IRAC 8.0 yum channel. The continuum and 
3.3 yum feature intensities are not correlated. 

We present model calculations which relate our measurements of the PAHs spectral energy distribution to the particles size 
and ionization state. Cation and neutral PAHs emission properties are inferred empirically from NGC7023 observations. PAHs 
caracteristics are best constrained in a line of sight towards the inner Galaxy, dominated by the Cold Neutral Medium phase : 
we find that the PAH cation fraction is about 50% and that their mean size is about 60 carbon atoms. A significant field-to-field 
dispersion in the PAH mean size, from 40 to 80 carbon atoms, is necessary to account for the observed variations in the 3.3 yum 
feature intensity relative to the IRAC 8.0 pm flux. However, one cannot be secure about the feature interpretation as long as the 
continuum origin remains unclear. The continuum and 3.3 yum feature emission process could be the same even if they do not 
share carriers. 

Key words. 

1. Introduction Red Array Camera (IRAC) on board Spitzer Space Telescope 

(SST), the sensitivity and angular resolution are available to 

Interstellar emission in the near infrared (NIR) traces the measure the NIR interstellar emission independently of stellar 

properties of the smallest dust particles known as polycyclic emission modeling, unlike with DIRBE. First determination of 

aromatic hydrocarbons (PAHs). Their presence in the dif- irac colors in Galactic fields was reported in the first round of 

fuse mterstellar medium (ISM) was inferred from photomet- Spj^^g^ publications by^ Bm^ from total power sky bright- 

ric measurements with the In^ra Red Astronomical Satellite ^^^^^^^ ^e undertake a more thorough study to quantify what 

(IRAS) (Bo ulangereta l. 1985), AROME dGiardet al. 1994) ^^^ ^^ jg^^^g^ ^^ p^H^ ^„^ tl^^ ^^^^^^^^ ^^^^ p^^icles with 

and the Diffuse Infra Red Background Experiment (DIRBE) jj^^ images of the diff^use Galactic emission. 
JDwek et alJll997B . First spectroscopic evidence was provided 
by th e Infra Red Telescope in Space (IRTS) (Tanakaetal. 



1996 1 and the Infrared Space Observatory (ISO) ( Mattila et al. Within PAH emission models, measurement of the 3.3 pm 



199a) for the Galactic plane. ISO succeeded to detect the feature is critical to constrain the PAH ionization state and size 



PAHs bands in spectra of high latitude ciiTus clouds for JLi & D l^aine"2001^. The existence of PAHs with a few tens 
/I > 5 jum JBoulanger et al. 2000) . Measurements of the of atoms was propos ed to account for the 3.3 pm emission 
shorter wavelength emission were still limited to bright ob- JLeger & Pugetlll984l) . Moreover, a continuum underlying the 
iects such as visua l reflection nebulae JVerstraete e t al."200lt 3.3 ^m fea ture has been detected in visual reflection nebulae 
Ivan Diedenhoven et alii2004) . For the first time, with the Infra JSellgren et al. 1983 ). and more recently in galaxies (JLu et alJ 

120031) ■ This continuum is not accounted for in PAH models and 



Send offprint requests to: its origin is still open : is it fluorescence emission from PAHs 

N. Flagey, e-mail: nicolas . f lageyOias . u-psud . f r or photoluminescence from larger grains ? 
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CVF spectrum {AjSA = 35 ^ 45) of the 



diffuse Galactic emission (for A^h — 10 cm" ) centered on the 
Galactic coordinates (26. 8, +0.8), completed for A < 5//m with 
the IRAC colors measured on the GLIMPSE field (diamonds). 
Horizontal lines represent IRAC filter widths. Dashed line : the 
gray body, with a color temperature Tc = 1105K and the 3.3 
fim feature, whose shape is extracted from the S WS spectrum 
of NG C7023 and intensity is the one measured bv lGiard et alJ 
1I1994I) . Vertical lines represents some eiTor bars on the spec- 
trum. 



In section 2, we present the IRAC data taken from the 
Galactic First Look Survey (GFLS) and Galactic Legacy 
Infrared Mid-Plane Survey Extraordinaire (GLIMPSE) as well 
as complementary ISOCAM/CVF spectra. We measure IRAC 
colors of the Galactic diffuse emission (section 3) and combine 
them with spectroscopic data in section 4 to separate the con- 
tributions of the 3.3 fim feature and the continuum to IRAC 
3.6 jum chann el. In sec t ion 5, we use an updated version of the 
model of De sert et alJ ( Il990l) - detailed in the appendix - to 
bring constraints on the PAH mean size and ionization state 
(section 6). We discuss the origin of the NIR continuum in sec- 
tion 7. 



2.5 
2.0 
1,5 
1.0 
0.5 
0.0 




1 12 




-2 



"^2.0 
-_ 1.5 
-_ 1.0 
^0,5 

_Jo.o 



Fig. 2. CVF spectra (A/6A = 35 ^ 45) of the diffuse Galactic 
emission centered on the Galactic coordinates (34.1,13.4), 
(299.7,-16.3) and (30.0,3.0), from bottom to top. Vertical Unes 
represents some errors bars on the spectra. 



We complement IRAC fields with ISOCAM/CVF spectro- 
scopic data covering the 5 to 16 j-im wavelength range at the 
positions listed in Table |4] over a 3' by 3' area. These obser- 
vations also point towards the diffuse Galactic medium, away 
from bright star forming regions. There is only one position 
common to IRAC and ISOCAM/CVF, centered on Galactic 
coordinates (26. 8, +0.8). We estimate the gas column density 
along this line of sight from th e H I Leiden/Dwingeloo sur- 
vey (Burton & Hartmann"l994') and the Columbia CO sur- 
vey (Cohen et al. 1986). The total column density is A^h - 
2.\QF^ cm"^ including 5.10^' H2 cm"^. For A'h, we assume that 
the emission is optically thin, and for A^h, we use the conver- 
sion factor 2.8 10^" H2 cm"^ per unit CO emission expressed 
in Kkms '. 



2. Observations 

2.1. Selected fields 

For this study, we use images from the GFLS and one field from 
GLIMPSE that span a range of Galactic longitudes and lati- 
tudes. These fields point towards the diffuse Galactic medium, 
away from bright star forming regions over path lengths which 
increase with decreasing Galactic longitudes and latitudes. The 
GLIMPSE field is centered on Galactic plane at a longitude / 
= 27.5° and extends over an area of 3° by 20'. It is a mosaic 
of 71x4 fields of 5' by 5' with an individual exposure time of 
2 seconds. The GFLS fields are centered on Galactic coordi- 
nates (l,b) = (254.4,+0), (105.6,+0.3), (105.6,+4), (105.6,+8), 
(105.6,+ 16) and (105.6,+32) and cover an area of 1° by 15'. 
They are mosaics of 3x12 fields of 5' by 5' with an individual 
exposure time of 12 seconds. 



2.2. ISOCAM/CVF spectroscopy 

The ISOCAM/CVF spectra are taken from the ISO Archive. 
We use the Highly Processed Data Products (HPDP) as de- 
scribed bv Boulanger et al. (2005) . The spectral resolution is 
between 35 and 45 for A between 5 and 16.5 fim. We pro- 
duce a mean spectrum by averaging all spectra over the 3' by 
3' field of view. In each spectrum the PAH features are visi- 
ble. The highest S/N spectrum, that centered on Galactic co- 
ordinates (26. 8, +0.8), is shown on Fig.^ The three others are 
plotted on Fig.|5] Error bars on the spectra, represented by ver- 
tical lines, were obtained by comparing spectra computed over 
distinct sub-areas. They are conservative estimates of the error 
bars as they may include true variations in the sky emission. 
The uncertainties are dominated by systematic effects (detector 
transients, zodiacal light subtraction) and are correlated over 
wavelengths (Boulanger et al. 2005i) . 
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Fig. 3. IRAC images (1° by 15') centered on Galactic coordi- 
nates (105. 6, +4) as viewed in 3.6,4.5, 5.8 and 8.0 //m channels 
(from left to right). Most of the point sources are masked (black 
dots on the images) and a median filter is applied. The dif- 
fuse Galactic emission is clearly visible in each channel. Color 
scales are different from one channel to another : the brightest 
areas ai-e in white (0.4, 0.3, 4 and 10 MJy/sr for 3.6, 4.5, 5.8 
and 8.0 yum channels), the faintest in black (0.1, 0.05, 1 and 3 
MJy/sr). 

2.3. IRAC images and processing 

All the IRAC images come from the Spitzer Archive ' (pipeline 
software version SI 1.0.2). We use the mosaiced images (post- 
bed). The diffuse emission is clearly visible in almost all IRAC 
data at a Galactic latitude below 16° (see Fig. but not at 
(105.6,H-32) which is located in a low column density region 
and that we use as an estimator of the noise. The surface 
brightness sensitivities are 0.0397, 0.0451, 0.154 and 0.165 
MJy/sr for IRAC 3.6, 4.5, 5.8 and 8.0 /im channels, as given 
by the Spitzer Sensitivity Performance Estimation Tool for a 
low background level ^. For the GLIMPSE field, the integra- 
tion time is shorter and the surface brightness sensitivities are 
0.409, 0.385, 1.03 and 0.886 MJy/sr for IRAC 3.6, 4.5, 5.8 and 
8.0 fim channels, for a high background level. 

We remove most of point sources by masking the pixels 
that are more than 3cr away from the image processed by a 
median filtering window of 21x21 pixels (1 pixel = 1.2"). Since 
we study the extended emission, we choose to apply a median 
filter (5x5 pixels) to the four IRAC channels data and keep 1 
pixel of every 3. This reduces the image size and increase their 
signal-to-noise ratio. The effective resolution of the images is 
then about 6 ". All images of a given field are projected on 



See http://ssc.spitzer.caltech.edu/archanaly/status/ 
See,http://ssc.spitzer.caltech.edu/tools/senspet/_ 



a common grid. We finally apply the photometric corrections 
given in Table 5.7 of the IRAC Data Handbook ^ for "infinite 
aperture" to all IRAC photometric results, since we focus on 
the diffuse emission extended over a significant fraction of the 
fields. After these corrections, the absolute calibration accuracy 
is about 5%. 



3. IRAC colors of the diffuse emission 

We measure the diffuse emission colors on the GLIMPSE and 
GELS fields by coiTelating the intensity in two IRAC channels. 
We discuss extinction correction and compute the ionized gas 
contribution to the colors. 



3.1. Correlations 

We measure the IRAC colors of the diffuse Galactic emission 
by correlating the brightness structure in each IRAC channel 
with the IRAC 8.0 ;um channel (see Eig.|5}. Some IRAC images 
present a strong intensity gradient along the long axis, that is 
obviously an artifact (see Eig.0}, which appears during the mo- 
saicing process and is due to bad dark-current correction. This 
gradient appears with different strengths on IRAC channels. It 
generally seems to be weak on IRAC 8.0 fim channel, rela- 
tive to the Galactic diffuse emission. On the other IRAC chan- 
nels, especially at 5.8 yum, the emission structure is sometimes 
dominated by a smooth gradient, that we take into account in 
the data correlation. We thus decompose each IRAC 3.6, 4.5 
and 5.8 fim images into a gradient, represented by a low order 
polynomial function of the long axis position, plus the emission 
structure of the IRAC 8.0 ;um channel. The uncertainties on the 
IRAC colors, associated with the gradient fitting, are estimated 
by looking at the variations of the color ratios with the order of 
the gradient. 

Eor each color, we iterate a linear regression, taking into ac- 
count the statistical noise per pixel as measured in the high lat- 
itude field (0.0268, 0.0222, 0.581 and 0.0639 MJy/sr for IRAC 
3.6, 4.5, 5.8 and 8.0 /im channel), and at each iteration, the pix- 
els that are more than 3cr away from the linear fit are masked 
for the next fit, where cr is the standard deviation of the dis- 
tance between the points and the straight line. Such pixels cor- 
respond to bright sources which are not removed by the median 
filter because they extend over an area comparable to the filter 
window. Eive iterations are sufficient to converge. 

3.2. Extinction correction 

Total hydrogen column densities are estimated from HI and 
CO observations (see section |5}. These column densities 
are converted into extinction in the IRAC channels combin - 
ing Av/Nh = 0.53 lO'^'cm^ fromlSayage & MathisI ( Il979l) . 
Ak/Av = 0.112 ('Rieke&Lebo fskvlll985h and the wave- 
length dependence of the extinction in IRAC channels de- 
termined from stellar measurements with GLIMPSE data 
llndebeto uw et alJ2005l) . We use their "Average" IRAC extinc- 
tions normaUzed to A^. The extinction corrections on the IRAC 



See_http://ssc.spitzer.caltech.edu/irac/dh/| 
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Fig. 4. Field centered on Galactic coordinates (105. 6, +8) as 
viewed by IRAC 5.8 /im channel before gradient correction 
(top), by IRAC 5.8 jjm channel after gradient correction (cen- 
ter), and by IRAC 8.0 //m channel before gradient correction 
(bottom). The small scale structures, visible on IRAC 8.0 fim 
channel before gradient correction exhibits on IRAC 5.8 fim 
channel only after gradient correction. 
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Fig. 5. Correlation plot of IRAC 4.5 yum channel versus IRAC 
8.0 fj.m channel for the GFLS field centered on Galactic coor- 
dinates (105.6,-1-4). The curves are iso-density contours (from 
the inside to the outside, the density of plot is at least 75%, 
50% and 25% of the maximum density). The straight line is 
the result of the linear fitting. The statistical noise per pixel, as 
measured in the field at (l,b) = (105.6,-1-32), is 0.0268, 0.0222, 
0.581 and 0.0639 MJy/sr for IRAC 3.6, 4.5, 5.8 and 8.0 fim 
channel. 



IRAC Aj channel, given the extinction r at both wavelengths 
and the observed color Ra^/a ,obs, measured in section lTTI : 



R 



Aj/Aj,int 



1 - eXp(-Ty) Ti 



(1) 



^3.6/8.0 is the only color to which we apply this correction. 
Differences in the extinction at 4.5, 5.8 and 8.0 //m are within 
eiTorbars (T4.5 = T5.8 = r^.o = 2.3 + 0.3). The 3.6 yum opacity is 
significant (T3 & = 3.1 at /? = 0°) but the /?3.6/8.o color coiTection 
is relatively small (1.24) because the opacities at both wave- 
length are close to each other. The extinction depends on the 
Galactic latitude. The b = 0° opacity gives the relevant extinc- 
tion correction since the GLIMPSE colors are dominated by 
low-latitude emission. In the data analysis we combine IRAC 
colors with AROME 3.3 //m feature observations and CVF 
spectroscopy. The extinction correction applied to the AROME 
measurements is consistent with ours (Giard et al. 1994) . The 
extinction correction for the CVF field, located at b = 0.8°, is 
negligible for A> 5 |/m. 

3.3. Ionized gas emission 

Free-free as well as gas lines emission might contribute to 
IRAC 3.6 and 4.5 fim channels, whatever the Galactic longi- 
tude. In order to obtain the contribution of the free-free emis- 
sion at 3.6 and 4.5 fim, we first measure this emission at 5 GHz. 
For the GLIMPSE field, we obtain the variatio n of the free-free 
emission from b - 1 ° to fe = 0° at 5 GHz from Alten hoff et alJ 
( 11979) . For the GFLS fields, we deduce t he free-free emis sion 
at radio frequencies from the Hq, emission ( lRevnoldsll992t) . We 
measure the H„ emission on the H,, Full Sky Map corrected by 
extinction bv Dickinson et al. ( 2003) . Then, we extrapolate the 
electrons emissio n (free-free and free -bound) from 5 GHz to 
NIR according to lBeckert et al.' {20Q(f). 

We finally add the contribution of gas lines. The fluxes 
from the main H recombination lines within the IRAC chan- 
nels (Pfs at 3.296jum, Pfy at 3.739 //m, Br„ at 4.051 yum, Pfg 
at 4.652/Ym, an d Pf,y at 7.46yum ) per un it Br^ emission are 
taken from the iHummer & StorevI ( 119871) Table for an elec- 
tron temperature and density of 7500 K and 10^ cm"^, and 
case B recombination (nebula optically thick to H ionizing pho- 
tons). At 3.6 //m, they represents 20% of the free-free emission, 
whereas at 4.5 fim, they double its contribution. The total con- 
tribution of the free-free and gas lines emission is thus about 
1% at 3.6 yum and 3% at 4.5 yum for most of GFLS fields, ex- 
cept the (105.6,-1-8) field, for which the contributions reach 3% 
and 11%. For the GL IMPSE field, the figures deduced from 
lAltenhoff et aP ( 1 1 979l) are about 7% and 12%. 



colors are not negligible for the GLIMPSE field, where the gas 
column density is the highest. For the other fields, the correc- 
tion is small compared to uncertainties. 

For the GLIMPS E field, our Ay is clos er than 10% to that 
given by flie maps of lSchlegel et alJ ( ll998l) wifli Ry = 3.l.We 
assume that the emitting dust is mixed with the grains responsi- 
ble for the extinction. The following relation gives the intrinsic 



3.4. IRAC colors 

Resulting IRAC colors are given in the first three columns of 
Table[T] For the GLIMPSE field, they are plotted, together with 
the CVF spectrum on Fig. [2 The ^5.8/8.0 ratios, given by the 
CVF spectrum (0.30) on the one hand, and by the GLIMPSE 
field (0.32) on the other hand, are in a good agreement. 



Table 1. IRAC color ratios. 



N. Flagey et al.; Origin of the NIR to MIR diffuse emission 



^3.6/8.0 



^4.5/8.0 



Vl.8/8.0 



V3.6/8.0,/sa( 



v3.6/8.0,fon( 



Nc 



SpitzerGFLS(105.6,+0.3) 


0.076 ± 


10 10-^ 


0.065 ± 8 10-3 


0.37 ± 5 10-- 


0.014 + 


2 10-3 " 


0.062 ± 7 10-3 


" 80 ± 20 


SpitzerGFLS(105.6,+4) 


0.059 ± 


8 lO-'' 


0.037 ± 5 10-3 


0.32 ± 4 10-' 


0.024 ± 


4 10-3 " 


0.035 ± 5 10-3 


" 56 ±12 


SpitzerGFLS(105.6,+8) 


0.094 ± 


15 10--^ 


0.050 ± 7 10-3 


0.26 ± 4 10-' 


0.047 ± 


10 10-3 a 


0.047 ± 7 10-3 


' 38 ±8 


SpitzerGFLS(105.6,+ 16) 


0.072 ± 


9 10-3 


0.046 + 5 10-3 


0.34 ± 4 10-2 


0.028 ± 


5 10-3 » 


0.044 ± 5 10-3 


- 52 ±12 


Spitzer GFLS (254.4,+0) 


0.068 ± 


8 10-3 


0.043 + 5 10-3 


0.37 ± 5 10-- 


0.027 ± 


5 10-3 " 


0.041 ± 5 10-3 


' 52 ±12 


Spitzer GLIMPSE (27.5) 


0.065 ± 


8 10-3 


0.052 ± 7 10-3 


0.32 ± 4 10-- 


0.016 ± 


2.10-3 b 


0.049 ± 6 10-3 


' 60 ±9 



Giardetal. (1994) 
Tanakaetal. (1996) 
Dweketal. (1997) 
Arendt et al. (1 998) 
Li & Draine (2001) 



0.016: 
0.016: 



2.10- 
4.10- 



0.0306 


0.0339 


0.0334 


0.0394 


0.0477 


0.0224 



0.278 



0.0264 ' 



0.0213 " 



The measurements uncertainties come from the linear fit- 
ting process, including gradient correction, and do not exceed 
3% on most of the IRAC colors. Extinction correction increase 
these uncertainties up to 10% whereas ionized gas corrections 
induce a negligible uncertainty. Taking into account the photo- 
metric accuracy (see section |3}, the final uncertainty is about 
13% on most of /?3.6/s.o, -^4.5/8.0 and /?5. 8/8.0, and reaches 16% 
for the GFLS (105.6,-h8) field. 

Our colors are averaged over large areas and do not give an 
account of the small scale dispersion across the IRAC fields. 
However, they already show strong variations from one field 
to another, especially /?3. 6/8.0 and /?4.5/8.o, whereas /?5. 8/8.0 does 
not vary that much around 0.3. 

4. 3.3 jum feature and continuum 

For the GLIMPSE field, AROME observations are combined 
with the IRAC colors to provide a spectrum of the diffuse emis- 
sion from 3 to 5 fim. We generalize this derivation to the other 
fields. 



4.1. Inner Galaxy spectrum 

Based on reflection nebulae observations, we interpret the 
GLIMPSE /?3.6/8.o and /?4.5/8.o colors with a PAH feature at 3.3 
yum and an underlying continuum. Spectroscopic observations 
of NGC7023 suggest that the feature and the continuum both 
contribute to the flux in IRAC 3.6 /im channel, whereas the 

flux in IRAC 4.5 jUm channel is dominated by the continuum 

J II — ^~l 

JSellgren et alJll983l) . The 3.3 yum feature has been spectro- 

scopi cally detected by IRTS in the inner Galaxy ( Tanaka et al. 
119961) . It ha s been measured p hot ometrically by t he AR OME 
experiment (^Giard et al. 1994) and B ernard et a l. (1994) show 
that the feature alone cannot account for the ISM emission in 
the DIRBE NIR channels. 

In reflection nebulae, the continuum is well described by 
a gray-body with a color temperature Tc - 1000 - 1500 K. 
This temperature is constrained thanks to photometric mea- 
surements on both sides of the 3.3 ;um feature. To interpret the 
GLIMPSE colors, we construct a NIR diffuse emission spec- 



trum with the spectral shape of the feature from the ISO/SWS 
spectrum of NGC7023 and a gray-body continuum. The fea- 
ture intensity is taken from the AROME measurement, taking 
into account the spectral response of their filters. We then fit 
the color temperature and intensity of the continuum to match 
the IRAC /?3.6/8.o and /?4.5/8.o colors. The resulting spectrum 
is shown in Fig. ^ The data does not constrain the spectral 
shape of the continuum. The combination of IRAC colors and 
the AROME feature only determine the ratio between the con- 
tinuum emission in the 3.6 and 4.5 fim IRAC channels. The 
color temperature of the gray body is Tc = 1 100 + 300 K, as 
computed from IRAC colors and AROME measurement. This 
value is in the range given by Sellaren et al. ( 1983) for reflec- 
tion nebulae, where the physical conditions are much different 
from the diffuse medium. The uncertainty on the color tem- 
perature is large because we only have two measurements to 
determine it. 

We can thus separate, for the inner Galaxy spectrum, the 
contributions of the 3.3 fim feature and the continuum to the 
IRAC 3.6 fj.m channel, which can be expressed as linear com- 
binations of the Ri.e/a.o and /?4.5/8.o colors (see Eq.|2land|3}- We 
define the contribution of the feature as the ratio between the 
flux of the 3.3 yum feature through the IRAC 3.6 /im channel 
and the flux of the PAH emission through the IRAC 8.0 jum 
channel, hereafter R3.6/s..Q,feat- 



^3.6/8.0,/«flr - ^3.6/8.0 - 0.95 X .R4.5/8.O 



(2) 



The contribution of the continuum to the IRAC 3.6 yum 
channel, is given by : 

Ri.6/&.0,cont = ^3.6/8.0 - ^3.6/8.0,/c«( 

= 0.95 X 7^4.5/8.0 (3) 

The corresponding values are listed in Table[n 

4.2. Additional fields 

For the GFLS fields, we cannot determine the value of the 
color temperature due to the lack of 3.3 fim feature measure- 
ments. Since the ratio between the continuum flux in IRAC 3.6 
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and 4.5 fim channels does not depend much on the color tem- 
perature, we assume that it does not vary among IRAC fields. 
Within this assumption, we use Eq.|2and|5] and give the cor- 
responding values of /^a.e/s.ojear and /?3.6/8.o,co,» for the IRAC 
fields in the two last columns of Table Q] It appears that the 
field to field variations of the 3.3 yum feature intensity are im- 
portant. There is a factor of 3.3 between the smallest values, at 
(/, b) = (105.6,-1-0.3) and (27.5,-t-O), and the highest value found 
for (/, b) - (105.6,-1-8). On the contrary, the field to field varia- 
tions of the continuum contribution /?3.6/8.o,conr are very weak. 
The mean value is about 0.04 and there is less than a factor 
of 1.8 between the extrema values. In all fields, the continuum 
has a strong contribution to the IRAC 3.6 jum channel. It ac- 
counts for one half of the flux in IRAC 3.6 /im channel at (I, b) 
= (105.6,-h8) and for more than 75% at (l,b) = (27.5,+0) and 
(105.6,H-0.3). There is no correlation between 3.3 jum feature 
and continuum colors. 

4.3. Comparison with DIRBE 

With the inner Galaxy spectrum (Fig. QJ, we c an compare 
our IRAC colors to thos e obtained with DIRBE (IDwek et alJ 
Il997t lArendt et al.lll998l) . We obtain these numbers, given in 
Table [0 by convolving the Galactic spectrum with the trans- 
mission curves of the instrument. For further studies, we give 
the conversion factors in Table |2 Our colors are significantly 
but slightly different from the DIRBE colors, especially for 
^3.6/8.0- However, this comparison does not lead to a unique 
conclusion. DIRBE colors corresponds to high-latitude emis- 
sion which may well differ from those measured on IRAC 
fields. The difference may also reflect systematic uncertainties 
in the DIRBE analysis associated with stellar subtraction. Last, 
but not least, our spectroscopic model may also contribute to a 
significant part of the difference. 

5. Spectral diagnostics of PAH size and ionization 
state 

The PAH emission s pectrum depends on both their charge state 
and size distribution. Li & Draine (2001) computed the diffuse 
ISM PAH spectrum with a model based on laboratory data. In 
Table \l\ we compare our colors to those expected from their 
model. Their model values are reasonably close to the colors 
which do not include significant continuum contribution. We 
have developed our own model, fo r several reasons. (1) The 
values listed by Li & Draine r2001) corresponds to a mixture 
between the Cold Neutral Medium (CNM), the Warm Neutral 
Medium (WNM) and the Warm Ionized Medium (WIM) in 



Table 2. Conversion factors from DIRBE or IRAS to IRAC 
fluxes, deduced from the Galactic spectrum. 



Instrument and channel 


Conversion factor 


DIRBE3.S/IRAC36 


0.917 


D1RBE4,)/1RAC4 5 


1.03 


DlRBEn/lRACso 


0.749 


IRAS 12/ IRACso 


0.863 
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Fig. 6. ^3.6/8.0,/eflf (solid lines) and ^7.7/11.3 (dashed lines) de- 
duced from the model as a function of the PAH mean size and 
G ^/f/n^. 



proportions 43%, 43% and 14% in mass, which are charac- 
teristic of high-latitude line of sight but do not apply to the 
low latitudes IRAC fields for which a significant fraction of the 
gas is molecular and thus are mostly sampling the CNM. The 
ionization state of PAHs along t hese l ines of sight might thus 
differ from that in Li & Draind ( 1200 ll) model. (2) As pointed 
out in previous sections, the 3.3 //m intensity measured through 
the R3.6/&.o,fear color Varies by a factor of 3 from field to field. 
We need to run a model over a grid of parameters to trans- 
late these variations in terms of PAH mean size and ioniza- 
tion state (see a ppend ix 1X1 for details). (3) In a recent study, 
tRapacioli et al.l ( 1200 5^ proposed a spectral decomposition of 
NGC7023 ISOCAM/CVF spectro-imaging data leading to dis- 
tinct emission spectra for cation and neutral PAHs. This opens 
the possibility to define the PAH properties from observations 
rather than l a borato ry measurements. Our model, an update of 
iDesert et alJ ( Il990l) . is based on this spectral decomposition, 
from which we derive the PAH cations and neutrals cross- 
sections, extending the work of Rapacioli et al. (2005). 

We run the model, coupled with a module that computes the 
PAH ionization fraction as a function of the PAH size. We apply 
this for various PAH size distributions and values of the ioniza- 
tion parameter G y/T/ne, where G is the integrated far ultravio- 
let (6-13.6 eV) radiation field expressed in units of the Habing 
radiation field, T is the gas temperatur e and »s is the electr onic 
density. The PAH size distribution of iDesert et al.l d 19901) has 
a mean size of 6 A or 45 carbon atoms, according to the rela- 
tion a - 0.9 V^ between the PAH size a in angstroms and the 
number of carbon atoms A^c- We vary the PAH mean size by 
changing the exponent of the PAH power law size distribution 
(standard value is -3) keeping fixed the values of the minimum 
and maximum PAH sizes (A^c = 20 - 180 or a = 4 - 12 A). 

We use two spectroscopic diagnostics to constrain the PAH 
properties : the ratio between the flux in the band at 7.7 and the 
band at 11.3 /vm (hereafter 7?7.7/n. 3) as well as Ri.e/a.ujeat- 
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5.1. R 



7.7/11.3 



as a tracer of PAH ionization state 



10000 



Using a Lorentzian decomposition that fits the PAH features 
between 6.2 and 12.7 jum, we calculate the ratio between the 
fluxes in the band at 7.7 and the band at 11.3 /vm (hereafter 
^7.7/11.3) for each PAH model spectrum. The central wave- 
lengths of the Lorentzian features are given as fixed inputs of 
the fitting process at 6.2, 7.6, 8.6, 11.3 and 12.7 fim, whereas 
the widths and amplitudes are set free. Fig. |6l shows ^7.7/11.3 
as a function of G 's/T/rif. and the PAH mean size. In the 
model and within the range of values we consider, this ratio 
depends much more on the PAH ionization than on the aver- 
age size : ne utral PAHs preserit a lower 7^77/11 3 ratio than ion- 
ized PAHs JDraine & Lill200lHBakes et alJl200lh . For a mean 
size of Nc - 45 carbon atoms, ^7.7/11.3 grows from 0.95 at 
G Vf /«e = 250 Ki/2 (.jn3 to 1 95 at G Vf /ne = 4000 K'^^ cm^ 
whereas it goes from 1 .0 to 2.2 for a mean size of A'c = 85 car- 
bon atoms, between the same values of GyfT/rie. This ratio 
traces the PAH ionization state. 

5.2. /?3.3/8.o,/ra? as a tracer of PAH mean size 

Since the model does not include the continuum underlying the 
3.3 fim feature, we immediately obtain ^3.6/8.0,/cfl/ by dividing 
the flux of the PAH model spectrum in the IRAC 3.6 jum chan- 
nel by its equivalent in the IRAC 8.0 /vm channel. Fig.|6lshows 
R3.6/s.o,fear as a function of the PAH mean size and G VT/ne. 
Unlike Rj. 7/113, /?3.6/8.o,/eaf depends on the PAH size as well 
as their ionization, and it is much more dependent on the PAH 
mean size when the PAHs are small. The 3.3 /urn feature be- 
come s fainter when the P AHs are big or ionized ( Draine & Li 
[20011 lBakesetalJl200lh . For a mean size of Nc = 45 car- 
bon atoms, Rj.e/s.Qjeat goes from 0.0175 when Gyff/rie - 
4000 K'''^ cm^ to 0.0375 when G Vf/we = 250 K'^^ ^.^^3 ^ and 
from 0.005 to 0.0125 for A^c = 85 between the same values 
of G y/r/rie. This ratio, when coupled with the previous one, 
constrains the PAH mean size. 

6. PAH size and ionization state across the diffuse 
ISM 

Combining ^7.7/11.3 and ^3.6/8.o,/«» (see Fig. |6j, we can now 
constrain the PAH mean size and ionization state for our in- 
ner Galactic spectrum, for which we have both these mea- 
surements. Adding the other IRAC colors /?3.6/8.o, ^4.5/8.0 and 
^5.8/8.0, we simultaneously constrain the shape and strength of 
the continuum. We also discuss the variations of the PAH prop- 
erties across the Galaxy. 

6.1. Inner Galactic diffuse medium 

For our inner Galactic spectrum, Ru/n.i - 1.5 ± 0.1 and 
Ri.b/&.o,feai - 0.016 ± 0.002. We compute a best-fit process 
over four parameters (A^c, G 's/T/rie, and the continuum inten- 
sity and color temperature) taking into account five constraints 
(the two ratios previously detailes, as well as the three IRAC 
colors R3.6/&.0, R4.5/&.0 and /?5. 8/8.0 given in TableQ. The result- 
ing x^ is plotted on Fig. 0and|8] The 3-cr error bars on the 
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Fig. 7. Iso-x^ contours deduced from a fitting process of 
^3.6/8.0, ^4.5/8.0, ^5.8/8.0, RiJiu.3 and Ri.b/s.Qjeai, plotted as a 
function of A^c and G Vr/«e. For each couple of A^c and 
G Vr/we, we take the lowest x^, depending on the intensity 
and color temperature of the continuum. Only the statistical 
dispersion on the colors is taken into account on this figure. 
Systematic and statistical dispersions account for the same in 
the resulting uncertainties on the parameters. 



parameters correspond to ;i'^ = 10. Fig. |5] shows the best-fit 
spectrum and corresponding physical parameters are listed in 
TableOl 

A PAH mean size of A^c = 60 car bon atoms i s large r 
by a factor of 1.3 than the mean size of IPesert et al.l lll99(l) . 
Along the inner Galaxy line of sight, half of the gas is molec- 
ular (see section |2j- A consistency check on the model is 
provided by the derived value of G y/T/rif. which should be 
close to t hat of the CNM . According to Li & Drains ( 200ll) 
and Wein gartner & Drain e (2001), in a radiation field G ^ 
1, G^/f/n^ = 250 - 380 K'^^ cm^ for the CNM, 2200 - 
2900 K'/2 cm^ for the WNM and 1000 K'^^ cm^ for the 
WIM. Considering that the radiation field increases from the 
inner Galaxy to the solar neighborhood by a factor of ~ 3 
(Sodroski et al. 1997), the ionization parameter for the CNM, 
near the molecular ring, is about 750 - 1 140 K' '^ cm^^, whereas 
the WIM and WNM are about 3000 K'^^ ^,^i ^nd 6600 - 
8700 K'^^ cm-'. Our Galactic spectrum is then in a good agree- 
ment with a CNM dominated medium, assuming that the mean 
CNM electron density and gas temperature are constant across 
the Galaxy. From the ionization parameter, we derive that the 
PAHs are half-neutral and half-cation. The contribution of the 
continuum in IRAC 3.6 yum channel, 70% ± 12% (see Fig.|8}, 
which corresponds to /?3.6/4.5,eonr = 0.90 ± 0.15, is in agreement 
with the one empirically determined in section 0] 

6.2. Field-to-field variations 

Sako n et alJ ( 120041) have measured 6.2, 7.7, 8.6 and 11.3 /vm 
features at different Galactic longitudes along the Galactic 
plane with IRTS spectroscopic data. They found that the 8.6 
and 11.3 yum features were systematically stronger relative to 
the 6.2 and 7.7 fim features in the outer Galaxy than in the 
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Fig. 8. Iso-x^ contours deduced from a fitting process of 
■^^3. 6/8.0, ^4.5/8.0, ^5.8/8.0, ^7.7/11.3 and /?3.6/8.o,/eaf as a function 
of the continuum contribution to IRAC 3.6 /vm channel and the 
PAH mean size. For each couple of A^c and continuum contri- 
bution, we take the lowest x^, depending on G Vr/ng ■ Only 
the statistical dispersion on the colors is taken into account on 
this figure. Systematic and statistical dispersions account for 
the same in the resulting uncertainties on the parameters. 
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Fig. 9. Bold dashed line : CVF spectrum of the diffuse Galactic 
emission (for A^h = 10^'cm"^) centered on the Galactic co- 
ordinates (26. 8, +0.8). Bold solid line : best fit resulting from 
the model, adding three components : PAH cations (dashed 
line), PAH neutrals [dash-dot line), and a continuum {dash-3- 
dot line). The parameters used for this fit are listed in Tabled 
The origin of the continuum is discussed in section0 

inner Galaxy. They suggest that PAH dehydrogenation or com- 
pactness may explain such band to band variations rather than 
ionization, even if they do not derive definite conclusions. Our 
CVF data allow us to address this question. 

In Table0] we list .R7.7/ 11.3 values for the diff'erent CVF lines 
of sight presented in section |2 The Lorentzian decomposition 
is made difficult on the CVF spectra by the weak S/N of some 
spectra and by the short range of covered wavelengths below 
6.2 and above 12.7 fim. In order to estimate the uncertainties 
of this measure, we try different methods, that differ by the 



Table 3. Best fit model output. Uncertainties result from statis- 
tical and systematic dispersions, which account for the same in 
the resulting uncertainties on the parameters. 



Model output 



Value 



Nc 

PAH+IPAH 

Continuum contribution to IRAC 3.6 



60±9 
1350±450 K'/2 cm^ 
42%±7% 
70%+ 12% 



Table 4. j^7 7/11 ^ from CVF lines of sight and from lSakon et alJ 



Line of sight 



Galactic coordinates 
(l,b) 



^7.7/11.3 



G26.8 
G34.1 
G30 
G299.7 



(26.8,0.8) 

(34.1,13.4) 

(30.0,3.0) 

(299.7,-16.3) 



1.5 
1.6 
1.3 

1.2 



tO.I 
tO.2 
tO.I 
tO.3 



Sakon et al. (2004) l=-8 


(-8,0) 


2.85±0.I7 


Sakon et al. (2004) 1=48 


(48,0) 


I.99±0.05 


Sakon etal. (2004)1=- 132 


(-132,0) 


I.42±0.40 


Sakon etal.(2004H= 172 


(172,0) 


I.65±0.37 



lorentzians parameters we fix or set free (position, amplitude 
or width). We give the ratios in Table |2 where we also list the 
ISakon et alJ ( 120041) values. Our two values for the inner galaxy 
(G26.8 and G30) are lower than those of Sakon et al. ( 2004) 
and there are no signs of a systematic Galactic gradient in the 
CVF data nor in Ri.b/s.ojeai color (see Table^. 

Within our model, the observed variations in the Rj.^/s.ojeat 
color of the Spitzer fields must trace variations in the PAH 
mean size. We refer to the solid lines plotted on Fig.|5] Since 
the stellar radiation field decreases from the inner to the outer 
galaxy, we consider that G y/T/ns^ is lower for the GFLS fields 
than for the GLIMPSE field (we assume that T and We are al- 
most co nstant from field to fiel d). Besides, we use the CNM 
value of iLi & Drain3 (I2OOII1 and lWeingartner & Drairi3 (l200lh 
G y/r/ne ^ 300 K'^^ cm^ as a lower limit because some of the 
gas must be in the WNM and WIM phases, where this ioniza- 
tion parameter is higher. This range of ionization factors allows 
us to constrain the mean size for each field, using the value 
of ^3.6/8.0,/fflr color. The derived PAH mean sizes are given in 
Table[n To estimate the error bars, we use /?3.6/4.5,co«/ - 0-9;^o2 
(resulting from the fit) in Eq.|2land|3] There is a factor of ~ 2 
between the two extreme values, at (105.6,-1-8) with 38 ± 8 car- 
bon atoms and at (105.6,-1-0.3) with 80 + 20 carbon atoms. The 
other PAH mean sizes are about 50 - 60 carbon atoms. What 
clearly appears is that we observe significant variations in the 
PAH mean size : the field at (105.6,-1-8) has much smaller PAHs 
than the field at (105.6,-1-0.3) and the GLIMPSE field. 

7. The origin of the NIR continuum 

The NIR dust continuum, first detected in visual reflec tion neb- 
ulae dSellgren et alJI 1983b and observed in galaxies JLu et alJ 
2003), is found to also exist in the diffuse ISM, where the ra- 
diation field is thousand times lower than in a reflection nebula 
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like NGC7023. This continuum accounts for 50% to 80% of 
the IRAC 3.6 /urn channel intensity in the diffuse medium and 
its field-to-field variations are weak, relative to the IRAC 8.0 
fiva flux (see Tableland Fig.Hll. We com p are our val ues to the 
one found for galaxies bv lLuetalJ ( l2003l) : ( 12004 . We esti- 
mate the contribution of free-free and gas lines emission to their 
IRAC 3.6 and 4.5 yum colors, using the Hq. to PAH emission 
ratio found by Roussel et alJ 



( 120011) . The|Lu| ( 12004 corrected 
colors R3.6/H.0 = 0.039 - 0.067 and R4.5/&.0 = 0.028 - 0.054 are 
thus in agreement with those we measure, even if they seem to 
be slightly smaller. 

The continuum observations raise two questions : what are 
the carriers ? what is the emission process ? 

It is not scattered light. With the scattering prope rties for 
dust in the diffuse ISM given bv lLi & Drainel (I2OOII) and the 
NIR interstellar radiation field derived from DIRBE NIR sky 
maps, we estimate the intensity of the scattered light per H nu- 
cleon. The scattered light accounts for 4% and 1 % of the mea- 
sured continuum in IRAC 3.6 and 4.5 /im channels in the inner 
Galactic spectrum. 

PAH fluorescence has been proposed to account for the 
continuum. We estimate the necessary photon conversion effi- 
ciency by dividing the number of photons emitted in the contin- 
uum between 2.5 and 5 jum by the number of photons absorbed 
in UV by PAHs. We assume that all of the UV energy absorbed 
by PAHs is re-emitted in the infrared and the mean energy of 
a UV photon absorbed by a PAH is 5.2ey ( Li & Draine 2001). 
We thus find t hat the necessary p hoton conversion efficiency 
is about 120%. lOordon et alJ (11^98) performed a similar cal- 
culation for the extended red emission (ERE) and obtained an 
efficiency lower-limit of 10% assuming that all the photons ab- 
sorbed by the dust are absorbed by the ERE producing material. 
Taking into account that PAHs are responsible for one quarter 
of the dust UV energy absorption, we find a NIR efficiency 
a factor of 3 higher than the ERE efficiency. If we consider 
a solid photoluminescence process - due to VSGs - instead 
of the molecular fluorescence, it leads to an equivalent photon 
conversion efficiency, because VSGs absorb almost the same 
energy as PAHs. In terms of energy, the molecular fluorescence 
or solid photoluminescence have an efficiency of about 8.5%. 

ISellgrenI (11984 have suggested that tri-dimensional grains 
of 45 to 100 carbon atoms that undergo stochastic heating like 
PAHs may be the carriers of the continuum in reflection neb- 
ulae. As mentioned in section \62\ the PAH mean size varies 
from A^c = 38 ± 8 carbon atoms to A^c = 80 ± 20 carbon atoms 
in the different IRAC fields we analyze. However, in our anal- 
ysis, the continuum is not correlated with the PAH mean size. 
[An & Sellgren (2003) found a systematic increase of the fea- 
ture to continuum ratio with increasing distance to the exciting 
star of NGC7023 . This could only reflect the dependence of the 
feature emission on the ionization state and not tell us anything 
about the continuum. 

The continuum emission questions the 3.3 fim feature in- 
terpretation. If there is an efficient fluorescence mechanism for 
the continuum, it could also account for the feature. Even if the 
carriers of the continuum and feature are distinct, we cannot 
discard the possibility that the emission process is the same. 



One cannot be secure about the feature interpretation as long 
as the continuum origin remains unclear. 

8. Conclusion 

We combine IRAC (GELS and GLIMPSE fields) and 
ISOCAM/CVF data to characterize the Near to Mid-IR 
Galactic diffuse ISM emission. Our results are as follows : 

- Extended diffuse emission is visible in most fields at the 
four IRAC wavelengths. We obtain IRAC colors by corre- 
lating the spatial distribution of the extended emission from 
one channel to another CVF spectroscopic data available 
on our main field, pointing towards the inner Galaxy, di- 
rectly demonstrates that the emission in the 5.8 and 8.0 //m 
channels comes from PAHs (mainly 6.2 and 7.7 jum fea- 
tures). On the same field, comparison with 3.3 fj.m feature 
photometric measurements show that PAHs accounts for 
only 25% of the IRAC 3.6 fim channel flux. A NIR contin- 
uum, previously seen in reflection nebulae, must be present 
to account for the remaining fraction of the IRAC 3.6 fim 
channel flux and the totaUty of the IRAC 4.5 yum channel 
flux. This decomposition is generalized to the other fields 
to obtain the 3.3 jum feature and continuum contributions to 
the IRAC 3.6 yum channel flux. 

- Among the fields we analyze, the PAH colors exhibit sig- 
nificant variations. On the contrary, the continuum remains 
almost constant, relative to the PAH flux in the IRAC 8.0 
yum channel and is in agreement with the values found for 
galaxies. 

- We interpret observed colors and their variations in terms 
of PAH mean size and ionization state. For this, we up- 
date the model of Desert et al. ( 1990) introducing different 
IR cross-sections for cation and neutral PAHs, determined 
from NGC7023 IRAC and ISOCAM/CVF observations. 
We use two spectroscopic diagnostics within the model and 
fit our main field observations with a set of 5 colors. The 
ratio between the 7.7 and the 11.3 jum features leads to a 
constraint on G y/T/rie which governs the equilibrium be- 
tween cationic and neutral PAHs. For the GLIMPSE field, 
we find G Vr/we = 1300 ± 200 K'^^ cm-' in agreement with 
theoretical expectations for a CNM dominated line of sight 
and the molecular ring radiation field. The cation fraction 
is about 50%. We also use the ratio between the 3.3 yum fea- 
ture seen by IRAC 3.6 fim channel and the PAH emission 
in IRAC 8.0 jum channel to derive constraints on the PAH 
mean size. Combining both ratios, we show that there is a 
significant dispersion in the PAH mean size, which varies 
among the IRAC fields, from 38 to 80 carbon atoms. 

- The continuum intensity is not correlated with the 3.3 //m 
feature, which implies that the continuum carriers might 
not be PAHs. It is not scattered light. A photon conversion 
efficiency of about 120% (energy conversion of 8.5%) is 
necessary to account for it by PAH fluorescence or VSG 
photoluminescence. 

The IRAC colors we have measured should be useful for 
further analysis of the diffuse emission in IRAC images. An in- 
teresting development of this work should be the analysis of the 
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IRAC colors and their variations on small spatial scales towards 
nearby clouds at high Galactic latitude to probe the evolution of 
PAH size distribution in interstellar clouds. If variations of PAH 
size distribution is confirmed, it could be correlated with varia- 
tions in the extinction curve and diffuse infrared bands (DIBs) 
in order to trace the contribution of small PAHs to these signa- 
tures. 

Appendix A: Determination of neutral and cation 
PAIH cross-sections 

In order to interpret IRAC colors and their variations in terms 
of PAH mea n size and ionization state, we update the model of 
IPesert et al.l (11990) to take into account the size depending ion- 
ization state of the PAHs. We thus introduce distinct emission 
properties for neutral and cationic PAHs. Immediately note that 
this model only gives an account of the PAH features emission 
but not of the continuum emission. The 3 to 13 //m absorp- 
tion cross sections of these two PAH forms are deduced from 
spectro-imaging CVF and IRAC observations of NGC7023. 

First, emission spectra /^(/l) and Iy{A) of pure PAH cations 
or neutrals are obtained through the following procedure. We 
start from the Rapacioli et al. ( 2005) spectra for A > 5/im based 
on their linear decomposition of the NGC7023 ISOCAM/C VF 
observations. With these spectra, we make a linear decompo- 
sition of CVF map which gives us weight-maps CAT{i, j) and 
NEUii, j) of the PAH neutrals and cations spectra contribu- 
tions (see Eg. lA. H and Fig. lA.H . CONST{i,j,A) accounts for 
emission of larger dust particles (VSGs). 



Iy(iJ,A) = ll{A)x NEUii, j) 

+ I^(A) X CAT{i, j) + CONS T{i, j. A) 



(A.l) 



With these two weight-maps, we decompose the difference 
map between IRAC 3.6 yum and IRAC 4.5 fim channel, accord- 
ing to Eq.|2l to obtain /+(/l) and Iy{A) at i = 3.6 yum. We also 
perform this decomposition on the IRAC 8.0 fira channel. 

The second step is to invert the emission spectra into ab- 
sorption cross sections of PAH neutrals and cat ions. We as- 
sume that the PAH size distribution is the one of IPesert et alJ 
(1199 0') model and we approximate the spectral distribution of 
NGC7023 radiation field by a black body of T,// = 17000 K 
(the effective temperature of the illuminating star). The inver- 
sion does not depend on field intensity. We iterate on the cross 
sections values until the model reproduces the emission spec- 
tra, and the intensity in IRAC 8.0 yum channel and 3.3 fim fea- 
ture to an accuracy better than 10%. 

The final amplitudes and widths for both PAH cations and 
neutrals are given in Table lA.ll The integrated cross-sections 
are given by the products of these two quantities. The cross- 
sections ratios for a given ionization state differ from those of 
Ilj & Draind (120011) by up to a factor of a few. The good cor- 
relation between PAH neutrals and the feature clearly appears 
on Fig. lA.H where contours are the 3.3 yum feature intensity, 
whereas PAH cations are closer to the star. The cross-sections 
of pure PAH cations and neutrals, for the Mathis radiation field 
and for the standard size distribution are shown in Fig. lA.2l 
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Fig. A.2. Neutral (solid lines) and cationic (dotted lines) PAH 
cross sections, for a 60 carbon atoms PAH, normalized to the 
7.7 yum value for neutral PAH. 

Table A.l. Amplitudes and widths for PAH cations and neu- 
trals that gives the integrated cross-sections, both normalized 
to the 3.3 yL( feature. 



A (Jim) 


^neutral 


Wneulral (cm-^) 


^colion 


Wcalion (cm-') 


3.3 


1.0 


20 


1.0 


20 


5.25 


0.013 


43 


0.82 


43 


5.7 


0.026 


43 


0.69 


43 


6.2 


0.401 


25 


24.0 


32 


6.85 


0.064 


65 


3.65 


85 


7.5 


0.129 


26 


9.60 


36 


7.6 


0.450 


31 


17.3 


26 


7.83 


0.206 


31 


14.4 


45 


8.6 


1.29 


20 


54.7 


25 


11.3 


3.34 


10 


28.8 


20 


12.0 


0.154 
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The model is coupled with a module that computes the 
cation, neutral and anion fractions for each PAH size from the 
local balance between the photo-ionization, the electron recom- 
bination, the electronic attachment and the photo-detachment 
processes as quantified by (Le Page et al. (2001)). The cation 
to neutral fraction depends on G yfT/ne, where G is the inte- 
grated far ultraviolet (6-13.6 eV) radiation field expressed in 
units of the Habing radiation field, T is the gas temperature 
and He is the electronic density. The fraction of PAH anions 
is computed to be small for diffuse ISM physical conditions 
(G Vr/ne ~ 100 - 3000). The small fraction of anions is as- 
sumed to have the same emission properties than neutrals. 
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Fig.A.l. NGC7023 weight maps of PAH cations {left) and neutrals (right), where the illuminating star is at the center and the 
black pixels code for the highest weights. The contours of Iv,33feature(i, j) are added in white (left) or black (right). Small scale 
contours (e.g. at (a, 6) = (21h01mn37s, 68°11'30")) are due to bright point sources, not taken into account in the fitting process. 
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